Thermal and ground-state entanglement in Heisenberg XX qubit rings 
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We study the entanglement of thermal and ground states in Heisernberg XX qubit rings with a 
magnetic field. A general result is found that for even-number rings pairwise entanglement between 
nearest-neighbor qubits is independent on both the sign of exchange interaction constants and the 
sign of magnetic fields. As an example we study the entanglement in the four-qubit model and find 
that the ground state of this model without magnetic fields is shown to be a four-body maximally 
entangled state measured by the A''-tangle. 



PACS numbers: Q3.65.Ud, 75.10.Jm 

Quantum entanglement is an important prediction of 
quantum mechanics and constitutes indeed a valuable re- 
source in quantum information processing Q|. Much ef- 
forts are devoted to the study and characterization of it. 
Very recently one kind of natural entanglement, the ther- 
mal entanglement |^-|^, is proposed and investigated. 
The investigations of this type of entanglement, ground- 
state entanglement [||^] in quantum spin models, and 
relations between quantum phase transition 

and entanglement provide a bridge between the quantum 
information theory and condensed matter physics. 

Consider a thermal equilibrium state in a canonical 
ensemble. In this situation the system state is described 
by the Gibb's density operator pT = exp{—H/kT) /Z, 
where Z =tr[exp (— iZ/fcT)] is the partition function, H 
the system Hamiltonian, k is Boltzmann's constant which 
we henceforth will take equal to 1, and T the tempera- 
ture. As pt represents a thermal state, the entanglement 
in the state is called thermal entanglement In a re- 
cent paper we showed that in the isotropic Heisenberg 
XXX model the thermal entanglement is completely de- 
termined by the partition function and directly related to 
the internal energy. In general the entanglement can not 
be determined only by the partition function. 

In this brief report, we consider a physical Heisenberg 
XX A^-qubit ring with a magnetic field and aim to obtain 
some general results about the thermal and ground-state 
entanglement. We consider a four-qubit model as an ex- 
ample and examine in detail the properties of entangle- 
ment. Both the pairwise and many-body entanglement 
are considered. 

In our model the qubits interact via the following 
Hamiltonian lllSi 
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i=l i=l 

where ai = {aix, cTiy, Giz) is the vector of Pauli matrices, 
J is the exchange constant and B is the magnetic field. 
The positive and negative J correspond to the antiferro- 
magnetic (AFM) and ferromagnetic (FM) case, respec- 
tively. We assume periodic boundary conditions, i.e., 
A^ + 1 = 1. Therefore the model has the symmetry of 



translational invariance. 

It is direct to check that the commutator [H, Sz] = 0, 
(rotation symmetry about the z-axis) which guarantees 
that reduced density matrix pi2 — tr3,4_...jv(/3T) of two 
nearest-neighbor qubits, say qubit 1 and 2, for the ther- 
mal state Pt has the form M 
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in the standard basis {|00), |01), |10), |11)}. Here \ij) = 
\i) (8) — 0,1), |0) (|1)) denotes the state of spin 

up (down), and Sz = X)fc=i '^jz/2 are the collective spin 
operators. The Hamiltonian is invariant under the trans- 
formation n^=i Sn^N-n+i, where Sij is the swap opera- 
tor for qubit i and j. This is reflection symmetry which 
implies that the matrix element z is a real number. For 
any operator A12 acting on qubit 1 and 2, we have the 
relation 



tri2(^12Pl2) = tri^2,---,A'(^12PT)- 



(3) 



Then the reduced density matrix pi2 is directly 
related to various correlation functions Ga[3 = 
(CTiaO-2/3) =tri,...,iv(o-iaO-2/3PT) [a = x,y,z) . Pre- 
cisely the matrix elements can be written in terms 
of the correlation functions and the magnetization 



M =iT:{Y.'i=i '^i^Pt) as 

=tri....,Ar(|00)(00|) 
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il + 2M + Gzz), 
(1-2A7 + C), 
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where M = M/N is the magnetization per site. In de- 
riving the above equation, we have used the translational 
invariance of the Hamiltonian. 

Due to the fact [H, Sz] = one has G^x = Gyy. Then, 
the concurrence [ |l4| quantifying the entanglement of two 
qubits is readily obtained as M 



1 



c 



0,|G, 



(1 + G, 



4A/2) 
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which is determined by the correlation function Gxx, Gzz, 
and the magnetization. In the Hterature there are a lot of 
results on correlation functions ]l5| in various quantum 
spin models, and they may be used for the calculations 
of entanglement. 

By using the translational invariance of the Hamilto- 
nian and the rotational symmetry about the z-axis we 
find an useful relation 



C = iU-BM)/{2J), 



(6) 



where U = U/N is the internal energy per site and U is 
the internal energy. From the partition function we can 
obtain the internal energy and the magnetization through 
the well-known relations 
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where (3 — 1/T. As a combination of Eqs.(^ and (|^) 
the correlation function Gxx is solely determined by the 
partition function. Therefore the concurrence can be ob- 
tained from the partition function and the correlation 
function Gzz- We see that the partition function itself is 
not sufficient for determining the entanglement. 

Except for the symmetries used in the above discus- 
sions there are also other symmetries in our model. Con- 
sider the operator A^^ = aix ® a2x ® ■ ■ ■ ® cfnx satisfing 
A^ = 1. We have the commutator [uiaOi+ia, h.x\ = 
and anticommutator [(7^2, Aj;]^ — 0. Then we immedi- 
ately obtain 

Gaa = tr{A^criaCr2a exp[-/3i7( J, B)]kx}/Z 
= tr{ai„a2a exp[-/3i7(J, 
M tr{A,ai, exp[-/?H( J, S)]A,}/Z 
= -tr{ai, exp[-/377( J, -B)]}/Z. 

These equations tell us that the correlation function Gaa 
and the square of the magnetization is invariant under 
the transformation B —B. From Eq.(H) and the in- 
variance of Gaa and we find 

Proposition 1. The concurrence is invariant under the 
transformation B ^ —B. 

The proposition shows that the pairwise thermal en- 
tanglement of the nearest-neighbor qubits is independent 
on the sign of the magnetic field. And it is valid for both 
even and odd number of qubits. 

In Ref. 1^, we observe an interesting result that the 
pairwise thermal entanglement for two-qubit XX model 
with a magnetic field is independent on the sign of the 
exchange constant J. Now we generalize this result to the 
case of arbitrary even-number qubits. 

For the case of even-number qubits we have 

H{-J,B) = AzH{J,B)kz, 

A^ = cri2 (g) (732 (g) • • • (g) CTJV-Iz- 



Th transformation A^ changes the sign of the exchange 
constant J. Note that the definition of A^ is different 
from that of A^;. It is straightforward to prove 

Gxx - tr(A,ai,(72.e-'^^('^'^)A,) 
= -tr{a,x<y2.e-f'"^-'-^\ 

Gzz = trCA.ai.as.e-^^^-^'^^A,) 
= tr(ai,a2.e-^«(-^^^)), 
M = tr(ai,e-'3^(-^^^)). 

From these equations we see that the absolute value 
|Ga;a;|, the correlation function Gzz, and the magnetiza- 
tion are invariant under the transformation J —J. 
Hence from Eq.(^ and the proposition 1 we arrive at 

Proposition 2. For the even-number XX model with 
a magnetic field the pairwise thermal entanglement of 
nearest-neighbor qubits is independent on the sign of ex- 
change constant J and the sign of magnetic field B. 

The proposition shows that the entanglement of AFM 
qubit rings is the same as that of FM rings. 

Now we consider the case of no magnetic fields. Then 
the magnetization will be zero and Eq.(^ reduces to 
G = i max [O, \U/J\ — Gzz — l] • We can prove that the 
internal energy is always negative. First, from the trace- 
less property of the Hamiltonian, it is immediate to check 
that 



lim U = lim 



-E„/T 



0, 



(8) 



where En is the eigenvalue of the Hamiltonian H. In this 
limit the correlation function Gxx and the magnetization 
M are also zero. Further from the fact that 

dU/dT = _ {Hf)/T^ = {^Hf/T^ > (9) 

we conclude that U is always negative. Then we arrive 
at 

Proposition 3. The concurrence of the nearest-neighbor 
qubits in the Heisenberg XX model without a magnetic 
field is given by 



„^ , 2max[0, -C//J- Gz^ - ] 
^ i max[0, U/J - Gzz - 1] 



for AFM, 
for FM. 



From the proposition we know that even for the case 
of no magnetic fields the partition function itself can not 
determine the entanglement. In the limit of T ^ 0, the 
above equation reduces to 



£;(0) 



^(0) 
(0) 



for AFM, 



i max[0, - 

imax[0,^-G^'i^-l] for FM 



(11) 



where E^^^ is the ground-state energy and G'f} is the 
correlation function on the ground state. Therefore the 
ground-state pairwise entanglement of the system is de- 
termined by both the ground-state energy and the corre- 
lation function gIz"* . Next we consider a four-qubit XX 
model as an example. 
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To study the pairwise entanglement we need to solve 
the eigenvalue problems of the four-qubit Hamiltonian. 
We work in the invariant subspace spanned by vec- 
tors of fixed number r of reversed spins. The subspace 
r = (r = 4) is trivially containing only one eigenstate 
|0000)(|1111)) with eigenvalue 4S(-4S). The subspace 
r = 1 is 4-dimensional. The corresponding eigenvectors 
and eigenvalues are given by 



1 4 

\k) = - ^ exp(— inA:7r/2)|n)o(fc = 0, . 



,3), 



(12) 



and 4Jcos(fc7r/2) + 2B, respectively. Here the 'num- 
ber state' \n)o = T"-i|l)o ,|l)o = |1000), and T is 
the cyclic right shift operator which commutes with the 
Hamiltonian H |16|. Due to the fact that H commutes 
Ax, \k)' = Ax\kjis also a eigenstate with eigenvalue 
4Jcos(fc7r/2) — 2B. The states \k) and \k)' are the so- 
called W states |l7| whose corresponding concurrence is 
1/2. We have diagonalized the Hamiltonian in the sub- 
spaces of r = 0, 1, 3, 4. Now we diagonalize the Hamilto- 
nian in the subspace r = 2 and use the following nota- 
tions: 

|l)i = |1100),|n)i=T"-i|l)i(,i = 1,2,3,4), 

|1)2 - llOlO), |m)2 = T"-i|l)2 (m = 1, 2). (13) 

The action of the Hamiltonian is then described by 

2 4 

H\n)i = 2J ^ \rn)2, H\m)2 ^2JY, I">i- (14) 

m—l n—1 

By diagonalizing the corresponding 6x6 matrix the 
eigenvalues and eigenstates are given by 

E± = ±4JV2, 1*^) -L \n)i ± ^ E l"^)^) ' 

Ei^O, = -i=(|1010)- 10101)), 

£^2 = 0, |*2> - ^(11100)- 10011)), (15) 

Es^O, 1*3) = ^(11001) -10110)), 

E4 ^ 0, 1*4) = i(|1100) -I- 10011) - |1001) - 10110)). 

From Eqs.([l^) and ( |l5|) all the eigenvalues are obtained 
as ±4B (1), ±2B (2), AJ±2B (1), -4 J± 2B (1), ±4 JV2 
(1),0 (4), where the numbers in the parenthesis denote 
the degeneracy. Then the partition function simply fol- 
lows as 

Z = 4 + 2 cosh(4V2/3 J) + 2 cosh(4/3B) 

+ cosh(4/3J)] cosh(2/3B). (16) 

From Eqs.(0), (|l|), and (|l|), we get the internal 
engery, magnetization, and correlation function Gzz , 



^ZU = 2jV2smh{AV2l3J) + 2Bsinh(4/3B) 

+2B[1 + cosh(4/3J)] sinh(2/3B) 

-h4Jsinh(4/3J)cosh(2/3B), (17) 
-ZM = 2 sinh(4/3S) + 2[1 + cosh(4/3J)] sinh(2/3B), 
ZGzz = 2cosh(4/3B) - cosh(4V2/3J) - 1. 

Then according to the relation (||) the correlation func- 
tion Gxx is obtained as 



GxxZ = -\/2sinh(4\/2/3J) 

- 2sinh(4/3J) cosh(2/3B). 



(18) 



The combination of Eqs.(|T^,(|T|), and (|) gives the exact 
expression of the concurrence. It is easy to see that the 
concurrence is independent on the sign of J and B, which 
is consistent with the general result given by proposition 
2 for even-number qubits. 

Figure 1 gives a three-dimensional plot of the concur- 
rence again the temperature and magnetic field. The 
exchange constant J is choosed to be 1 in the follow- 
ing. We observe a threshold temperature Tc — 2.36338 
after which the entanglement disappears. It is interest- 
ing to see that the threshold temperature is independent 
on the magnetic field B for our four-qubit model. For 
two-qubit XX model the threshold temperature is also 
independent on _B |^. We further observe that near the 
zero temperature there exists a dip when the magnetic 
field increases from zero. The dip is due to the energy 
level crossing at point B^ = 2{V2- 1) = 0.82843. When 
B increases form B = 2 near the zero temperature the 
entanglement disappears quickly since there is another 
level crossing point Bc2 = 2 after which the ground state 
becomes |1111). We also see that it is possible to increase 
entanglement by increasing the temperature in the range 
of magnetic field B > 2. 

Now we discuss the ground-state entanglement (T = 
0). When B < Bd, the ground state is l^"-) with the 
eigenvalue i?*^"-' = — 4V2J. It is direct to check that 
ciz^ = —1/2. Then according to Eq.(pl]), the concur- 
rence is obtained as C = V2/2 - 1/4 = 0.45711. When 
Bel < B < Bc2 the ground state is |fc = 2) and the corre- 
sponding concurrence is 1/2. When B > Bc2 the ground 
state becomes |1111), so there exists no entanglement. 

Thus far we have discussed the pairwise entanglement 
in both the ground state and the thermal state. Now 
we discuss the many-body entanglement in the ground 
state. Recently Coffman et al ||T^ used concurrence to 
examine three-qubit systems, and introduced the concept 
of the 3-tangle, as a way to quantify the amount of 3- 
way entanglement in three-qubit systems. Later Wong 
and Christensen generalize 3-tangle to even-number 
A''-tangle. The TV-tangle is defined as 



Tl,2,...,Af = KV'lo'ly ® <T2y 



CTNy\r)\\ (19) 



with |i/)) a multiqubit pure state. The A^-tangle works 
only for even number of qubits. 
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For the ground state \^-) we have = \^-), and 

aiy (g) a2y (E)<T3y (E) aiy\-9 -) = 1"$ -) , (20) 

i.e., the ground state is an cigcnstatc of the operator 
cij/(8'(T2j,(8'(73j,(E)cr4y . Therefore we find Ti^2,.. .,4 = 1, which 
means that the ground state has four-body maximal en- 
tanglement. For another two ground states when varying 
the magnetic field it is easy to check that ti.2,...,4 = 0, 
which means that the ground state have no genuine four- 
body entanglement. Note that the ground state \k — 2)' 
has pairwise entanglement but no four-body entangle- 
ment. 

In conclusion we have found that the pairwise thermal 
entanglement of nearest-neighbor qubits is independent 
on the sign of exchange constants and the sign of mag- 
netic fields in the XX even-number qubit ring with a 
magnetic field. For determining the concurrence we need 
to know not only the partition function but also one cor- 
relation function Gzz- For the four-qubit model we ob- 
serve that there exists a threshold temperature which is 
independent on the magnetic field. The effects of level 
crossing on the thermal entanglement and ground-state 
entanglement are also discussed. Finally we find that the 
ground state is a four-body maximally entangled state ac- 
cording to the potential many-body entanglement mea- 
sure, A''-tangle. 
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FIG. 1. The concurrence against the temperature and mag 
netic field. The parameter J — 1. 
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